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We describe hadronization events, using the AdS/CFT Correspondence, which display many of the 
quahtative features expected in QCD. In particular we study the motion of strings with separating 
end points in a back-reacted hard wall geometry. The solutions show the development of a linear 
QCD-like string. The end points oscillate in the absence of string breaking. We introduce string 
breaking by hand and evolve the new state forward in time to observe the separation of two string 
segments. A kink associated with this breaking evolves to the end points of the string inducing rho 
meson production. We explicitly compute the rho meson production at the end point. 



INTRODUCTION 

The AdS/CFT Correspondence [H 121 H [4] provides a 
new tool for the study of strongly coupled gauge theories 
and has been extended to cover both gauge and quark 
[SmiTllHlElITOlIIIllig degrees of freedom. Hadroniza- 
tion, the conversion of quark pair production into jets 
of hadrons, is traditionally a thorny problem in QCD. 
It is therefore interesting to model this process with the 
Correspondence. 

The combination of asymptotic freedom and confinement 
in QCD ensures that the road from quark production 
to jets, as in a high-energy accelerator such as LEP, 
is broken up into three overlapping stages. This is the 
statement of factorization. For a recent review, see [13]. 
The first, pair production, typically occurs at high en- 
ergies and so the relevant degrees of freedom are a bare 
quark/anti-quark pair. As they fly apart, these quarks 
begin to radiate soft gluons and quark/anti-quark pairs 
- the process known as showering - so that the relevant 
degrees of freedom are dressed quarks. Showering can be 
computed at high energies in QCD with effective theories 
including soft collinear effective theory. After the dressed 
quark/anti-quark pair have propagated a distance of or- 
der 1/Aqcd apart, the strong coupling and confining dy- 
namics of QCD become relevant: the quark pair forms a 
highly excited colourless hadron that can and does frag- 
ment. This last conversion is known as hadronization. 

The middle process of showering has been studied in dif- 
ferent contexts via the Correspondence. A number of 
authors [14, il5j have computed the evolution after the 
injection of R charge into the AT = 4 Super Yang-Mills 
theory. The A/" = 4 theory is strongly coupled at all 
scales and so there is no suppression of the emission of 
large transverse momentum - rather than seeing jets, the 
events fill the entire space surrounding the initial pair 
creation. 



In the AdS/CFT correspondence, strings with their ends 
tied to a D7 brane represent a quark anti-quark pair in 
the dual field theory. The dual quarks have a constituent 
mass that goes to zero in the limit that the D7 brane fills 
the entire AdS space. The evolution of such pairs in the 
massless quark limit with separating end points was first 
studied in [16 . They worked in the quenched M = 2 
gauge theory obtained by placing probe D7 branes in 
AdSs X . Both the endpoints and the string between 
them fall indefinitely into the radial direction of the AdS 
space representing the spreading of baryon charge and 
energy respectively in the dual field theory as the state 
evolves. These states are anisotropic and their evolu- 
tion naturally describes the showering of flavour. Hadron 
(rho meson) production can only occur on the D7 world- 
volume through the motion of the endpoints of the string. 
The production for these states was computed [16 and 
jet-like structures were seen. 

Much attention has focused on heavy quark propagation 
at finite temperature [I9l |20] . The A/" = 4 theory is per- 
haps more like QCD in its non-conformal high tempera- 
ture phase which is described by an AdS-Schwarzschild 
geometry [3 . These computations have also been in- 
spired by attempts to reproduce jet quenching observed 
in heavy ion collision data at RHIC. 

In this paper we wish to return to study the second and 
third phases of hadronization described above at zero 
temperature. To find a dual gravity description that be- 
haves like QCD one must seek a non-conformal version 
of the AdS/CFT Correspondence that at least incorpo- 
rates confinement. There has been a considerable body of 
work in this direction. We will work in the context of de- 
formed AdS geometries [2ll[2a[23|2l[23[2i[271|2i[^ 
which have the benefit of being 3+1 dimensional at all 
energy scales and having naive perturbative dimensions 
for a class of operators in the UV (of course the UV is 
not perturbative but conformal, strongly coupled, and 
highly supersymmetric). Simple deformations such as 



the introduction of masses for the A/' = 4 matter content 
that break the supersymmetry to A/" = 2 [22], AT = 1, 
[231, [H US] or less [11 [271 |M1 have been studied. 
These duahties are not as weh understood as the original 
Correspondence but a number of features seem generic. 
In particular, the geometries dual to confining theories 
develop a wall structure in the interior that stops fields 
from penetrating to energy scales (radial distances in 
AdS) below the mass gap of the theory. Such a repulsive 
wall leads to linear potentials between quarks since they 
are described by the ends of a fundamental string in the 
geometry - typically the centre of these strings fall on to 
and lie close to the repulsive wall with their energy then 
just scaling with their length. 

In fact, the imposition of a simple hard wall cut off 
on AdS is now frequently used to model a mass gap 
[301 ISB [32] . That will not be sufficient for us since we 
wish to evolve the strings onto the wall and therefore need 
the repulsion to be described by the geometry. Given 
that these dualities are not fully understood we will take 
the simplest example of a deformation with such a hard 
wah. So called "dilaton ffow" geometries [271 [281 [29] have 
a non- trivial dilaton profile in a deformed AdS space. 
The dilaton, or on the field theory side the gauge cou- 
pling, blows up at some infra-red scale generating the 
wall (effectively the scale Aqcd)- The wall is repulsive 
to strings and also to the D7 branes they are attached 
to, which dynamically induces both confinement ^3 [33j 
and chiral symmetry breaking [TOl [H]. Such a geome- 
try is clearly a decent first stab at describing a QCD-like 
gauge theory. 

Having found a model for such a theory, we follow 
the program of \16^ and evolve strings in this geome- 
try. In particular, we consider numerical solutions of the 
Polyakov string action for strings with ends tied to a D7 
brane and point-like initial conditions to represent the 
pair creation of a quark anti-quark pair. We allow the end 
points to separate quickly leaving a string falling into the 
geometry between. It is fairly easy (though computer- 
time consuming) to numerically follow this evolution with 
pre-built numerical integrator functions such as ND Solve 
in Mathematica. As one might expect, in the hard wall 
setup the string falls onto and bounces off the wall in the 
radial direction - the centre of the string then oscillates, 
falling on to the wall repeatedly. Whilst this happens, the 
endpoints move apart and a longer string grows, taking 
up the kinetic energy near the end points. Formally, we 
study this process at infinite N and so there is no string 
breaking; instead, the string continues to evolve. Since 
the string can only extend so far, the endpoints eventu- 
ally stop and their motion reverses. The string bounces 
in the three dimensional space. The model therefore cor- 
rectly reproduces expectations from field theory in the 
large N limit. 



For real QCD string breaking is presumably greatly 
favoured, occurring essentially as soon as there is of or- 
der Aqcd worth of potential energy in the string. In 
our gravity description, string breaking is captured by 
stringy 1/N effects: there will be some transition ampli- 
tude between our initial string and a final state and the 
actual evolution will involve a sum over all possible inter- 
mediate states, including those with two string segments. 
Attempts have been made to compute this breaking prob- 
ability in [34, 35 , 36 , 37 which we do not add to here. We 
choose instead to study the evolution after string break- 
ing, so we simply break our strings by hand and consider 
particular states with two string segments. In particular, 
we insert a static quark/anti-quark pair and then evolve 
these new particular strings. We simulate the inserted 
pair with a static string that falls straight down into the 
AdS space from the D7 brane. At a time of our choosing 
(usually once the initial string grows to a size of order 
1/Aqcd), we break the initial string in the middle by 
hand and join each half to a static segment. We then 
follow this broken string as it evolves. 

This configuration again evolves as one might expect - 
the static endpoint (anti-quark) is dragged by the fast 
moving endpoint (quark) that it has been attached to 
and the endpoints (quark/ ant i- quark pair) of the orig- 
inal string are free to separate to infinity. The broken 
strings have a kink at the point where the static string 
was attached which evolves to the endpoints. The kink 
then causes a rapid jerk in the motion of the inserted 
static endpoint (anti-quark). 

In the AdS/CFT Correspondence, mesons made of the 
quarks are associated with open strings tied to the D7 
brane |9j - in fact, in the large N and strong coupling 
limit, rho mesons are the dominant such mode, described 
by a dual gauge field on the D7 world volume. The string 
endpoint is a source for that gauge field. The rapid jerk of 
an endpoint provides a mechanism for copious production 
of the rho mesons, in the same way that an accelerating 
charge radiates in classical electromagnetism. We show 
how to compute this production. 

The rho meson and each of its radially excited states is 
associated with a normalizable mode in the radial direc- 
tion on the D7 brane. The full production of the gauge 
field on the D7 world-volume can therefore be computed 
mode by mode. In fact the radial dynamics simply en- 
codes the mass of the physical state and the strength 
with which a particular source couples to that mode. For 
a single mode, the calculation is just of the usual form 
for a moving source coupled to a massive vector field in 
4d and is straightforwardly computed. A static quark 
(string endpoint) has a cloud of mesons around it with 
their density decaying with distance as an exponential of 
the mass. When an endpoint accelerates, it also radiates 
waves which correspond to the production of rho mesons 



in the hadronization event. We explore this production 
for the rho meson and its excited states. 

The solutions we present for moving string solutions pro- 
vide a qualitative understanding of how hadronization 
may occur in QCD or related theories. We do not in 
this paper wish to make any claims that these computa- 
tions are numerically accurate for QCD, but we hope that 
the mechanisms revealed will provide thought for the fu- 
ture modelling of hadronization in QCD. We discuss the 
prospects in our final section. 
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Figure 1: The basic D3-D7 system showing the coordinate 
labelling we use - note the D3 lies at r = 0. We also show 
two possible moving string configurations - one lies off the 
plane of closest approach between the D3 and the D7 and 
lies in a 3d space, the string with end points separating along 
that plane lies in just 2d. We study moving strings of the 
latter type. 

GEOMETRIES AND QUARKS 

The vacuum structure of AT = 4 SU {N) super Yang-Mills 
gauge theory is encoded in the AdSs x geometry 



(1) 



where = AnggNa'^^ Qs is the string coupling, and 
a' sets the string length. The dilaton is constant and 
there is a four form C(4) = /R^dx^ A ..dx^. Note that, 
in these coordinates, r (the radial direction in the six 
dimensional vq space) corresponds to the gauge theory 
energy scale (r = is the infra-red). The AdS space 
has a boundary at r = oo and the dual gauge theory 
"lives" on that boundary with a metric conformal to the 
Minkowski space metric dx^j^^. 

To introduce a single quark fiavour (really a single M = 2 
hypermultiplet), we place a D7 brane lying in the x^ — x^ 
and — directions (see Figure 1). In pure AdS, such 
a brane lies fiat with rl^r^ = r^. The separation ro sets 
the scale of the constituent quark mass (m^ = ro/27ra'). 



related to the energy of an open string with one end at the 
bottom of the D7 brane and the other at the AdS horizon 
r = 0. In contrast, a fundamental string with each end 
on the D7 brane represents a quark/anti-quark pair and 
the interaction energy between them. As shown in [TT], 
for very heavy static quarks (i.e. with their ends tied 
to a D7 brane at r = oo), the dual string hangs further 
and further into the bulk as the quarks are separated. 
The total energy of the state scales as one over the quark 
separation. 

To provide a gravity description of a confining gauge the- 
ory, we will consider a geometry that is asymptotically 
AdS (the UV contains the same degrees of freedom as the 
AT = 4 theory) but which has a back-reacted hard wall in 
the IR. The simplest such example we know is a dilaton 
controlled deformation [271 IM] • The metric in Einstein 
frame for this deformed geometry is 



ds" = ^A\r)dxl^^ 



(2) 



with 



3/2 



and the four form is as in the pure AdS solution. 



(3) 



Note that the dilaton (dual to the gauge coupling) and 
metric have a singularity at r = r^^ which can be loosely 
interpreted as the scale at which the gauge coupling di- 
verges, Aqcd- In [38j it was argued that the best inter- 
pretation of this geometry might be as dual to the low 
energy theory below some UV cut off where some highly 
irrelevant supersymmetry breaking but SO (6) preserving 
coupling becomes important in the A/" = 4 gauge theory 
(a scalar mass term is an example of such an operator). 
For our purposes in this paper though, the geometry is 
simply playing the role of a back reacted hard wall to 
include confinement and the precise physics is unimpor- 
tant. 




Figure 2: A sketch of a D7 embedding in the back reacted 
hard wall geometry for a quark with zero current mass. The 
bending of the D7 induces a constituent mass (non-zero 
minimum D3-D7 separation) for the quarks. 



Quarks can generically be added to the theory by in- 
cluding a D7 brane as in the AdS case. The hard wah is 
repulsive to the D7 brane and so induces chiral symmetry 
breaking through the formation of a quark condensate - 
the details are discussed in [lOl |TT] . In Figure 2 we sketch 
the form of the D7 embedding in the geometry for the 
case of zero bare quark mass - the D7 is bent so that 
there is never a zero length string between the D3 and 
D7: the quarks have acquired a constituent mass. 

As discussed in |^ a string with ends tied to a D7 
brane, corresponding to a quark pair, is repelled by the 
hard wall at r = and the string lies close to that wall 
- there is a linear potential between the quarks. 



Also, we label the timelike world sheet coordinate as r 
and the spacelike as a (i.e. dPcj = drda). 



The equations of motion are 



daX^dbX 



N d 



(6) 



PN 



0. 



Note that the derivative in the final term acts on r]^^ 
depending on the choice of S. 

The open string boundary condition applied at the string 
end points a = <j* is 



N 



(7) 



Conventions 



Throughout this work, we consider fields living on four 
different spaces: the full lOd geometry, the worldvolume 
of a D7 brane, the worldsheet of a fundamental string, 
and on the Minkowski spacetime of the dual field theory. 
To simplify some of this chaos, we choose some conven- 
tions for how we label indices and these fields. First, 
we will use the upper-case latin letters (M, A^, ..) to re- 
fer to lOd spacetime indices (e.g. GMNdx^ dx^ is the 
metric on the lOd space), while lower-case latin letters 
{a^b^ ..) refer to worldsheet and worldvolume indices. Fi- 
nally, we use lower-case Greek letters (/i, i^, ..) to index 
the Minkwoski directions in both the lOd geometry and 
D7 worldvolume (i.e. x^, /i = 0, 1, 2, 3). 



However, if the string is attached to a D-brane localized 
at = c^, then this condition is replaced in the direc- 
tions transverse to the brane by the Dirichlet condition 
X 



M 



There are also world sheet constraints from variation with 
respect to rjah are 

daX^dkX'^GMN = Ivab^'GMNdcX'^ddX^, (8) 

which, with the form of riab above become 

X-X' = 0, (9) 



X^ + T?X'^ = 0, 



(10) 



where X indicates drX and X' denotes d^X. 



STRINGS IN PURE ADS 



A study of string evolution in pure AdS was presented in 
[16]. Here we repeat that framework, which we will use 
below, and summarize the form of the results. 

We work with the Polyakov action in Einstein frame 



1 



Alia' 



d ere 



0/2^ 



r]r]''^daX''''dbX'^G 



-AT. 



MAT, 



(4) 



where rjab is the world sheet metric, labels the world- 
sheet coordinates, X^ are the embedding functions, and 
Gmn is the background metric. One is free to use the 
symmetries of the theory to fix the world sheet metric 
and we will pick the form 



Vab 



-S(x,r) 






1 



(5) 



We pick and choose the form of the stretching function S 
problem by problem in order to keep our numerics stable. 



We will also need to compute the energy of the strings 
that we generate. That energy can be determined from 
the lOd space-time stress energy tensor 



5Sf 



-G 5Gmn{,x) 



^rj'^'daX'^dbX 



(11) 

The conserved energy is then given by the integral over 
the 9d space of Tq , 

^"d^ J d^'^e'f/^y^v"''9aX°dbXo5it-X°). (12) 



If we rewrite the delta function as 

5(t - n(a)) 



5{t-X°) 



\drX0{T4a),a)y 



(13) 



where r^{(j) is the solution to X^(r*(cr), a) = t, we can 
then integrate over r to find 



27ra' 



\drX°\ 



T = T:t, (cr) • 



(14) 



Separating Quark Solutions 



In [T^ falling strings in AdS with endpoints separating 
in the Minkowski directions were studied - those so- 
lutions describe states with massless quarks in the basic 
M = 2 theory described by the D3-D7 system. For mass- 
less quarks, the D3 lies within the D7 and the strings 
with ends on the D7 are free to come arbitrarily close 
to the D3s at r = 0. The solutions they found had the 
endpoints separating in x and, at late times, falling in 
r along geodesies. The centre of the string fell faster in 
r leaving a curved string between the endpoints. They 
showed that kinks in the initial condition of the string 
were quickly smoothed out by the inflation-like growth 
of the string. 

We now wish to begin to make contact to QCD where, 
although the current quark masses are small (for the 
light quarks), they nevertheless obtain a large constituent 
mass from the chiral symmetry breaking dynamics of the 
theory. For this reason, we will work with the D3-D7 
system in a configuration where the D3 and D7 do not 
coincide (i.e. the dual quarks are massive). Later, we 
will consider the setup where we embed the D7 in the 
deformed AdS of Eq. 

As shown in Figure 1, the motion of a string ending on 
the D7 is generically a three dimensional problem. The 
string can move in the Minkowski directions x^ and the 
r^~^ directions on the D7 whilst the centre of the string 
may droop into the remaining two r^~^ directions of the 
geometry. For simplicity, we will restrict the motion to 
two dimensions - that is, we will look at strings lying at 
the point of closest approach between the D3 and the 
D7, with ends separating in the x^ directions. There 
is no hard computational block to performing the more 
generic computation, but we believe we can pick out the 
crucial physics within the simplified set-up. 

So that our analysis can be compared to that in p!6], 
we will follow their example and define an inverse radial 
coordinate 

p4 

' (15) 



The D7 brane embeddings are more complicated in these 
coordinates, but we will only work with strings at the 
point of closest approach to the D3, zq. The strings we 

) + r? + r? =0. The con- 



-. The 



study therefore lie at + + + r\ 
stituent mass of the quarks is then — ^ / • 

^ y 2t^ oi' zq 

AdS boundary is now located at z = and the IR (or 
the D3 branes) is at z = oo. The string then falls in 
the 3d space parameterized by the Minkowski directions 



t,x = X and radial coordinate z with metric 



The general initial conditions for such a string are pa- 
rameterized by specifying the embedding functions 
and z and their derivatives at world sheet time r = 
as a function of a. As an example, consider the initial 
condition 



t(cr,0) = 0, x(cr,0)=0, z{a,0) = zo. 



(17) 



Note that this initial condition corresponds to a configu- 
ration where the string is initially contracted to a point. 
Since X'(cr, 0) = 0, the first world sheet constraint Eq. 
91) is satisfied and the second Eq. (10) is simply = 0. 



Of course, the point-like initial conditions Eq. (17) are 



not sufficient to specify the subsequent string evolution. 
We must also specify the r-derivatives at r = consis- 
tent with the constraint equations Eq. (|8|. Consistent 
point-like initial conditions are then specified by two free 
functions, x(cr, 0) and i((j, 0). The last derivative, t(cr, 0), 
is fixed by the constraints to be 



(18) 



As an example, we consider strings with initial deriva- 
tives 

i(cr, 0) = Czo sin cr, x(cr, 0) = I^zq cosa, (19) 

where the string endpoints are located at a* = 0, tt. This 
set of initial conditions has two free parameters C and D 
that essentially set the string's momentum in the z and 
X directions respectively (the total string's x momentum 
obviously vanishes - D sets the size of the momentum in 
the X direction of either half of the string). Strings gener- 
ated from large D and small C will have end points that 
quickly separate in the x direction, so that the dual field 
theory state contains a back-to-back quark/anti-quark 
pair with large momenta. 

The string generated from these initial conditions has 



some energy. Using Eq. (14), and the initial value of 
t(a,0). 



sin^ a ^ D'^ cos^ a, (20) 
we find that the total energy of this string is 



sin^ a -\- cos^ a. 



(21) 



(16) 



where we have normalised the stretching function S by 
S(^o) = 1- The constant C lets us choose the initial 
speed of the centre of the string in the z direction whilst 
D (in the large D/C limit) picks the total energy of the 
string. 

With these initial conditions we can solve the three equa- 
tions of motion (IgI) provided we also impose the end point 



boundary conditions 

t'(cr*,r) = 0, x'(cr*,r) = 0, z{a^,T) = zq. (22) 

Note that our choices of initial time derivatives above 
are also consistent with these boundary conditions. Fi- 
nally, since the initial conditions satisfy the constraint 
equations Eq. ([8|, the evolved string satisfies them at all 
times. 

A judicious choice of the stretching function allows ND- 
Solve in Mathematica to follow the evolution for a con- 
siderable time. Here we pick 



2 / 2 



(23) 



To check the consistency of the solutions along the time 
evolution we monitor the total energy of the configuration 
using Eq. ( 14 ) - for all our solutions energy is conserved 
at, at least, better than the 1% level through the evolu- 
tion. For the simple solutions in AdS the conservation is 
much better than 1%. 

In Figure 3 we show the evolution of such a configura- 
tion with C and D = 1 and zq = 0.5R. As the end points 
separate, the string droops down into the bulk of AdS. 
Increasing C only serves to push the middle of the string 
down into AdS faster. Overall, this picture is consistent 
with the usual expectations from static strings in AdS 
and with the results described at zero mass in [1^. In 
particular, the fact that the string can penetrate indef- 
initely into the bulk of AdS refiects the conformal 1/r 
potential between two quarks. Of course, this behaviour 
is closer to the asymptotically free regime of QCD than 
the confining phase. We now turn to modelling the latter. 



STRINGS WITH A HARD WALL 



We will now study the motion of strings in the deformed- 
AdS geometry Eq. Q to represent quark/anti-quark pair 
production in a confining gauge theory. As for the strings 
we studied in the last section, we will again set the string 
in motion in the x direction whilst localised in the r^~^ 
directions at the point of closest approach of the D7 to 
the D3 brane - see Figure 2. The embedding function of 
the D7 brane will not therefore play any role in the com- 
putation other than determining that closest approach 
point, Zq. In practice, even zq would be a phenomenolog- 
ical parameter that would need to be fitted to the cur- 
rent quark mass - for the qualitative analysis here we set 
^w/zo = 2 but the precise value is not important (if one 
made the ratio very large then the quarks would become 
heavy relative to the scale Aqcd ^ ^I^w representing a 
large bare and constituent quark mass). 





Figure 3: The evolution of a string with end points on a D7 
brane at = 0.5i? in pure AdS with the initial conditions 
Eq. ( 17|19 ) and parameters C — D — \. The top plot 
shows the string worldsheet evolution and the bottom a 
series constant time slice shots of the strings motion at 
t = 0.4i?, 0.8i?, 1.2i?, and 1.6i?. 

We can now consider the evolution of strings in this ge- 
ometry using initial conditions like those we used above 
in pure AdS - that is, we look at a point-like string with 
separating end points. For simplicity, we will take almost 
exactly the same initial conditions. We maintain the con- 



ditions in Eqs. (17), (19) but to satisfy the constraint Eq. 



(10) with the new metric, we need a new initial condition 



for t(cr,0): 



t(cr, 0) = (C7A2(zo)) sin^ cr + L)2 ^^^2 (24) 

Since A{z) rapidly becomes unity away from the wall 
z = Zyj^ these strings are generated from essentially the 
same initial conditions as in the pure AdS analysis. We 
use the same stretching function S as in the AdS case 



here, Eq. (23). We plot the resulting string evolution in 
Figure 4. 



singularity 



qq creation 



D7- 

Figure 5: Sketch of the insertion of a quark anti-quark pair 
into the string evolution to represent string breaking. 



STRING BREAKING 




Figure 4: String evolution in a hardwall geometry. Here 
we set C — 2.5, D — 0.5 in the initial conditions Eqs. 
( 17|19 ), and choose zq — 0.5^1^ to set the quark mass. The 
top plot shows the world sheet evolution. The bottom 
plot shows time slices through the world sheet at times 
t — 0.2zw, 0.5zw, O.Szw, and l.lzw 

The solution plotted in Figure 4 shows many of the prop- 
erties one would expect. In particular, as the centre of 
the string approaches the hard wall, it is repelled and a 
straight section of string builds against the wall - this is 
the formation of the naive QCD string. We are work- 
ing at infinite N where string breaking is forbidden and 
so, rather than breaking, the string continues to evolve. 
The centre of the string begins to bounce off the hard 
wall and starts oscillating between the position of the D7 
brane and the hard wall. Meanwhile, as the potential 
energy in the string grows and the quarks separate, the 
quarks' kinetic energy begins to be sucked into the string 
- they slow. Eventually the quarks are brought to a halt 
and the string between them begins to contract, reversing 
the quarks' motion until they pass through each other. 
This oscillation will continue indefinitely in the absence 
of string breaking. 

We conclude that the inclusion of a hard wall does indeed 
begin to move the description of quark pair production 
closer to the expected behaviour in QCD. In the next 
section we will discuss how to include string breaking 
and thereby allow two separated jets to emerge. 



Our analysis of hadronization so far lacks one crucial in- 
gredient with respect to QCD - there is no string breaking 
and therefore no hadronization! String breaking is an in- 
herently 1/N effect and so absent from the AdS/CFT 
Correspondence in the tractable limit we study. Some 
progress has been made in studying this process though 
- if the string lies within the D7 worldvolume where it 
can break directly the rate is computable [34[ |35] |36l [37] 
and is of order 1/N. When the string lies away from the 
D7 there must be a quantum fluctuation that brings the 
string back to the D7 where it can break - estimates for 
this process can be found in [36l [37] . Here we do not wish 
to add to these computations. Our instinct instead is to 
assume that in QCD the string will break with essen- 
tially probability one as soon as there is of order Aqcd 
of potential energy stored in the string. We can model 
this on the gravity side by breaking our string at the 
corresponding time. 

We can insert a quark/ ant i- quark pair into the evolution 
of our strings by hand. An unpaired quark corresponds to 
a string stretched from the D7 brane to the D3 branes. A 
static solution of this form with a new choice of stretching 
function (S = A{z)/A{zo)) is 



t(r,a) 



TvAizo) ' 



(25) 



z(r, cr) = ^ {zo - Zs)a^ 2zs 



where x is constant. This solution satisfies the equations 
of motion, the constraint equations and the boundary 
conditions for all r and a. It is a straight string stretching 
from zq (at <j = tt) through the point Zg (at <j = 7r/2). 

We can therefore split our string in the centre and com- 
plete it to the D7 brane with a straight string segment as 
shown in Figure 5. As an example we will split the string 
solution shown in Figure 4 at the point where the centre 
of the string starts to bounce off the wall {t = 1.13z^) - 
the QCD string has roughly just formed at this point in 
the time evolution. At the splitting time ts (r^ in world 
sheet time) when t(r5,7r/2) = tg and z(rs,7r/2) = Zs we 
use the new initial conditions for the range 7r/2 < <j < tt. 



Strictly inserting the extra string length does not con- 
serve energy for the half string piece - our solutions below 
do have sufficient initial energy that the straight string 
is slightly less than a 10% correction to the total energy. 
We don't expect the non-conservation to have any great 
effect on the qualitative behaviour of the solutions we dis- 
play. In fact since we only study the evolution of one half 
of the initial string configuration one could imagine that 
some asymmetry in the distribution of energy between 
the two halves might be induced in the string breaking 
in any case. One could also splice out a straight string 
section in the middle of the initial string before inserting 
the vertical string pieces back to the D7 brane but we 
again would expect to see little qualitative change in the 
behaviour of the solutions. 

The time evolution of the half string segment can now 
be followed and we show such a numerical solution in 
Figure 6. The solution again behaves in accordance with 
naive expectations - the fast moving end point of the 
string continues to move, whilst two kinks in the string 
induced by the breaking propagate to the end points. 
When the kink arrives at the static end point, it is jerked 
(along with the entire string segment) in the direction of 
motion of the fast end point. Similar evolution of broken 
strings may be found in Following the evolution 

without further string breaking leads to oscillations of the 
end points in the segment's centre of mass frame, which 
separates infinitely far away from the other broken string 
segment. We also plot in Figure 7 the motion of the two 
end points of the string through this evolution - the string 
end points are special because they source a gauge field 
on the D7 brane world volume which corresponds to rho 
meson production. We will compute this production in 
the next section. 
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Figure 7: Plot of the string end point motion for the 



configuration in Figure 6. 




0.0 0.5 1.0 1.5 2.0 

X 



Figure 6: We show a string world sheet evolution after 
imposing string breaking. The initial condition involves half 
of the string from Figure 4 at time t = l.lSzw broken and 
extended by a straight string back to the D7 brane at a = 7r/2 
to cr = TT. The top plot is again the world sheet evolution 
and the bottom plot slices taken at t = lAzw, l.Szw^'^-Szw, 
and 3.2zw 



RHO MESON PRODUCTION 



Rho mesons 



To understand how the string solutions above radiate 
energy into hadronic modes, one must study the elec- 
tromagnetic theory on the surface of the D7 brane. For 
simplicity, we will assume that the D7 of Figure 2 has 
been pushed out sufficiently far away from the hard wall 
in the geometry Eq. Q that we may treat the lOd geom- 
etry as AdS. Further, we assume that the D7 is essentially 
fiat. That is, we will approximate the D7 embedding with 
the flat embedding of a D7 in pure AdS down to r = ro 
with induced metric 



P[G]^g=P^dx^ 



i?2 



i?2 



A^^aT^^dp'' + p''dnl). (26) 



Here, p is the radial direction on the world volume of 



the D7 so that p^ 



r\. Although these 



approximations may seem a little crude we will see that 
the p dependence of the problem enters essentially just 
through the mass of the mesonic states which can there- 
fore be easily altered to any more complicated set up. 

The end points of the string act as electrically charged 
sources for the gauge field that lives on the world volume 
of the D7 brane. The equation of motion for that gauge 
field follows from the variation of the electromagnetic 
action on the brane, 



i j d^x^FatF-'^ I cfx^gfAa. (27) 



The variation of the gauge field gives both the equation 
of motion and the boundary action. 



EM 



1 



(28) 

SO that the equations of motion are just Maxwell's equa- 
tions. 



1 



J 



and there is a boundary action (at p = 1/e ^ oo) 

SSk^y = - [ d^xdnsV^SAaF'"'. 
Jp=i 



(29) 



(30) 



Let us ffist understand the solutions in the absence of a 
source. 



The propagating modes of the gauge field arrange them- 
selves into multiplets of the 50(4) isometry group of the 
S^. The resulting Kaluza-Klein fields each map to dif- 
ferent operators of the gauge theory; in particular, the 
singlet on the maps to a conserved baryon current [9] . 
To study this current, we therefore give the solutions of 



Eq. (29) without sources for the modes that only have 



non-zero and are singlets on the S . We impose the 
gauge choice V^A^ = 0. The equation of motion in the 
absence of sources is then [9] 



-^Mp'dpA,) - (^2 f ,2)2 V(4)-4m = 0. (31) 



where V^4^ = rf^^d^dy is the scalar Laplacian in 
Minkowski space and P is a second-order differential op- 
erator. Fourier transforming in the Minkowski directions, 
we can write the equation above as an eigenvalue equa- 
tion in the radial coordinate with eigenfunctions fn{p) 
given by 

f {p)-I 2Fi(-n,-n-l,2; V/rg) 

where the In are normalization constants, n = 0,1,2..., 
and eigenvalues 



M2=4(n + l)(n + 2)-^ 



(33) 



e^fn{p)e'^ 



The solutions to Eq. ([31| are therefore given by the 
modes A., 



with kz 



These states have a discrete mass spectrum and are iden- 
tified with the rho mesons of the dual gauge theory. The 
factor of Tq /i?^ ~ / A indicates that the meson masses 
are much smaller than the quark mass at large 't Hooft 
coupling. Moreover, the fn (with appropriate choices of 
the In normalizations) are orthonormal functions (sub- 
ject to the weight factor w) 



dp w /„/, 



W 



(34) 



{p^+riy 

The corresponding choice of normalization is given by 

ll = 2{n + l)(n + 2)(3 + 2n). (35) 

We are using the solutions appropriate for the M = 2 D3- 
D7 configuration here rather than those for a bent brane 
as in Figure 2. However, as we will see, the holographic 
directions only enter into our final radiation computation 
through the masses they endow the four dimensional rho 
mesons and the value of J^, the normalization of the wave 
functions. In the more complicated case, one could sim- 
ply switch the spectrum and normalizations as appropri- 
ate. 



Green's Functions 



Boundary Data 



To observe the emission of rho mesons by the string end 



points we wih solve (29) by means of a Green's function 



for the field A^. For the minimum-energy configuration 
quark (a string connecting the D3 and D7 at their point 
of closest approach) the endpoint lies at p = where the 
volume of the three-sphere is zero and hence the Green's 
function is a constant on the three-sphere (the equiva- 
lent of the fn functions for R-charged states fall to zero 
at p = 0). This implies that there is no production of 
i?-charged rho mesons associated with non-trivial spher- 
ical harmonics on the S^. For a more generic string mo- 
tion such states would be produced. Moreover, since the 
sources do not move in the r^~^ four-plane, both the 
radial and angular components of the source current 
vanish and thus Ap and Ai (the components of the gauge 
field along the S^) also vanish in this gauge. We can 
therefore consider only the Minkowski components of the 
Green's function. 

Having chosen the Lorentz gauge, that Green's function 
satisfies 



The near-boundary behaviour of the gauge field is related 
to the one-point function of the dual conserved baryon 
current in the field theory. In particular, that one-point 
function is given as 



(J^(x^)) = lim 



UGRA 



oM^(x-,l/e)' 



(39) 



where 5'sugra is the on-shell bulk gravity action and the 
bulk gauge field A^ is the singlet mode on the S^. Using 
the variation of the bulk action in Eq. ( 28 ), the boundary 
current is simply 

{.^ixn) = -limyrj'^''d,A,ix-',p)\,=y,. (40) 

We can therefore write a bulk-to-boundary Green's func- 
tion that relates the bulk source to the boundary current. 
In particular, we write 



(fx'^g^'^,{x^;x')f {x'), (41) 



1 



r.3 M 



(36) 



where we define the bulk-to-boundary Green's function 
Q as 



where V is the differential operator defined in Eq. (31). 



Since the equation of motion for the gauge field in the 
presence of our source is given by 



(37) 



the full solution for an arbitrary current distribution 
follows from the convolution integral 



g^^,{x-;x'\p') ^ ^2(-l)"/„/„(p')G^,^,(a;^x"'), 

n 

(42) 

where the fn are the eigenfunctions in Eq. ( 32 ) and Gn 



is the 4d Green's function for a massive vector as before. 
The factor of 2(— 1)"^/^ comes from the insertion of the 
near-boundary expansion of the /^'s. 



f„(p) = (-l)("+l)/„4+o(/,)-^ 



(43) 



A^{x) = j d^x'^G%{x,x')ii^,,,f'{x'). (38) 



into the form of the boundary current in Eq. (40) 



The actual current distribution will be localized on the 
worldline of the string endpoint and will take the form 
= qx^5^{x) where the dot represents differentiation 
with respect to proper time. 

In order to obtain the Green's function, let us ex- 
pand in the basis of eigenfunctions describing the rho 
mesons used in Eq. (32) so that (p, x^; p', x'^) = 
Y.nfn{p)fn{p')G^p{x'',x"'). Inserting this form into 
Eq.(^6|, multiplying by fm and integrating over all 



space we find that the four-dimensional functions Gn 
are just the Green's functions for massive vectors in 
Minkowski spacetime with masses corresponding to the 
rho meson masses. 



Retarded Potential 



We have now reduced the problem to solving for each 
mode Gn the retarded potential for a massive field in fiat 
space. The retarded potential takes the form [4^ 



1 

47r 



9{t-t'){5{a) + V{cT)e{-a))5^^.. (44) 



Here we use the Synge world- function a = \iripy{x — 
x'Y{x — x'Y . The non-singular part of the solution is 
given by V{cf) = —-^^Ji[Mn^/—2a) where Ji is the 
Bessel function of order 1. 



static String End Point 



Radiation From String End Points 



As a first example of using this formalism we will com- 
pute the baryon density around a static quark. Consider 
such a charge at x = and at p = (r = ro), the point of 
closest approach on the D7 brane. We will concentrate 
on the temporal component of the gauge field which 
is dual to the operator '^7^'^, the quark density. 

One seeks to evaluate the integral over the past trajectory 
of a point source moving with a constant speed in a 'static 
gauge' given by x' = /3t' . Doing the spatial integral using 
the fact that the source is located at a point in the space- 
like dimensions leaves one with the integral 



47r 



/ dt' 9{t-t') {5{a) + V{a)e{-a))^ 



dr d£_ 
dr 



(45) 



where we have used the fact that fn{p = 0) = In- 

The time component of the four- velocity is actually can- 
celled by a Lorentz factor coming from the splitting of 
Minkowski spacetime into space-like sections when in- 
tegrating along the particle worldline. The two con- 
tributions to the integral are easily computed (for the 
non-singular piece due to the massive field using the in- 
tegration variable u = Mn\/—2o-)) giving the correctly 
Lorentz-covariant expression (the 7 is the usual boost 
factor) 



47r y^72(a;-/?t)2 + y2 + ^2 



(46) 



In the rest frame of the point source this reduces to the 
usual Yukawa form. The full solution is a sum over modes 
weighted by the fn normalizations see Eq. (35). 



There is a rapid rise in these normalizing factors with 
n which is due to the end point of the string being a 
delta function. Away from N ^ 00 one would expect the 
string to have some width and the expansion to truncate 
at some intermediate n. In any case this rise is not faster 
than the exponential fall off of the solutions so the physics 
away from the source is still dominated by the lightest 
modes. The Green's function converges for all |x| > due 
to the exponential factor in the Yukawa potential of each 
partial wave. The behaviour is dominated by the lighter 
modes at distances comparable to the Compton wave- 
length of the lightest mode. We interpret this Green's 
function as the 'dressing' of an isolated quark by a cloud 
of mesons. Holography gives the relative amounts of each 
of the excited states in the cloud. 



We now have a framework in which the emission of 
mesons can be modelled using the techniques of classical 
relativistic wave equations. The retarded Green's func- 
tion is straightforwardly integrated over the past world- 
line of an accelerating endpoint to give eg. for a particle 
moving in the x-direction (the u variable is as defined in 
the preceding section) 



(J0(X^))„ 



2(-l)"J^.9 
47r 



t-t' (cr=0) -^{X-X' (cr=0)) 



Jiju) 



t-t' 



- {x—x') 



(47) 



We will again plot the baryon number density which is 
holographically encoded by the sum over the J^. It may 
be noted that the plots we obtain give the superposition 
of radiated baryon density and the static baryon density 
associated with the probe quark. An elementary pre- 
scription is available for computing the reaction force on 
the probe quark due to the radiation (by differencing the 
advanced and retarded potentials) but this is not what 
we are interested in here (it involves a negative count- 
ing of the non-causal advanced potential and so would 
not produce a plot resembling meson emission). In holo- 
graphic scenarios (large A^) the force exerted on the quark 
by the dynamics of the colour fiux tube far exceeds the 
reaction force from meson emission anyway. In the case 
of an instantaneous acceleration it is possible to subtract 
the appropriate static and boosted solutions inside and 
outside of the particle's light cone but we do not apply 
this here. 



Massless Meson Limit 



In the strict A ^ 00 limit the meson masses are very 



small relative to the string mass (see ( 33 ) . At least for the 



lightest members of the tower, it is therefore interesting 
to compute the radiation into a massless gauge field on 
the D7. For this case cr = \ {-{t - t'f + (x - x'{t')f). 
The first term in (47) then gives 



47r 



^ 2{-irilq r da5{a) 

47r J (^(^t'-t) + {x'-x)^) y^^) 

2{-iril 1 

47r t-t'{cr=Q)-Xo{x-x'{cr=Q)) 

This is straightforward to evaluate for acceleration kicks 
such as those we found for the string end points in Figure 



7 above. For example the static end point is accelerated 
quickly to a constant speed. As an example form for 
the function x'{t') that describes a stationary particle 
accelerating to a final speed a we take 



x'{t') 



ab 

TT 



-at' 



tan 



(49) 



b here controls the time interval over which the accelera- 
tion occurs - we plot some sample trajectories in Figure 




i.o 

(491) used 



to 



parameterize the motion of an accelerating point source. The 
parameter a controls the final speed and is set to a = 0.2 
here, b controls the time scale of the acceleration and the 
plots show b = 0.8 (top), b = 0.3 (middle) and b = 0.05 
(bottom). 




Figure 9: The radiation of light mesons by a quark given an 
impulse in the positive x-direction, shown in the z = plane. 
The plot shows the density of the emitted (radiative part of 
field) and bound (boosted static part) mesons. The top plot 
is for a terminal velocity of 0.2c and below is 0.6c (in both 
plots the parameter b — 0.2). 



It is a simple matter to plot the resulting wave induced. 
Emission is typically a spherical shell radiating from the 
point of acceleration - there is an SO (2) symmetry in the 
7/, z coordinates so we shall plot the intensity of the wave 
in the x, y plane at 2; = 0. Examples of the gauge field 
produced are shown in Figure 9. The radiative piece is 
visible along with the 'hill' of the boosted static poten- 
tial. A clear, narrow emission wave is observable. For 
larger values of the final speed a the forward emission 
is typically enhanced relative to the backwards emission, 
and the overall emission is greater. For smaller acceler- 
ation times (smaller b) the wave front simply becomes 
narrower. For the accelerations of the string end points 
in Figure 7 we expect precisely such emission of the lower 
mass members of the mesonic tower. 



Massive Meson Limit 



For members of the meson tower with masses close to 
the quark mass (very high n at large 'tHooft coupling) 
we must compute the non-singular term in (47) which 



involves numerical integration of the Bessel function. In 
Figures 10 and 11 we show the effect of increased meson 
mass on the radiated mesons (the meson mass should be 
compared to the inverse time over which the string end 
point is accelerated). As the meson mass is increased we 
find emission of the more massive states are suppressed. 

In the massive case the waves are dispersive and produce 
an interesting pattern which is not just a wave localized 
on the light-front. The 'wavy' emission of meson density 



can be considered to be arising from quantum-mechanical 
interference effects. 

In principle we could sum over the emission of all of the 
meson states. At large 'tHooft coupling though there are 
many states lighter than the quark mass so the result 
would be unilluminating. The precise form of the meson 
masses and the coefficients are also model dependent. 
However, we believe that the computations we have made 
show how in principle the radiation could be computed 
and give a good understanding of the generic features of 
that meson radiation. 




Figure 10: The radiation of massive mesons by a quark given 
an impulse in the positive x-direction, plotted along the x- 
axis. The plot shows the density of the emitted mesons. The 
meson masses increase through the plots as 0,10~^, 10~^, |, 1. 
Note the background static field peak becomes narrower as 
the mass is increased. 
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Figure 11: Emission of massive vector mesons (m = |)- The 
parameters a — b — 0.2. 



LESSONS FOR QCD 

We have explored how a hadronization event happens in 
a gauge theory that has the degrees of freedom of A/" = 4 
super Yang-Mills in the UV but a deformation that leads 
to a back reacted hard wall and confinement in the IR. 
This theory, of course, is not QCD but the generic picture 
that emerges may have some lessons for the construction 
of a phenomenological model of hadronization in QCD. 

In particular we have suggested a picture in which the 
initial quark anti-quark pair separate, growing a string 
between them that dips into a holographic radial direc- 
tion. Initially the string's energy represents a 1/r poten- 
tial between the quarks reminiscent of the asymptotically 
free regime of QCD. In actuality the A/" = 4 theory is 
strongly coupled - one should ignore any radiation from 
this string (such as open string (glueball) radiation in 
the bulk) that would be present in the AdS/CFT Cor- 
respondence since that would be the wrong physics for 
QCD. The string then encounters a hard wall in the ge- 
ometry and spreads out along it forming a traditional 
QCD-like string with energy growing with its length. We 
have proposed that this string will break quickly once 
there is sufficient energy in it to pair create a quark pair, 
breaking the string. The two sub-strings then separate, 
radiating energy through their end points as rho mesons. 
We expect the majority of the hadronization particle pro- 
duction to come from this radiation rather than repeated 
string breakings. 

It is interesting to compare this picture to that of the 
Lund string model [41 which is one of the leading de- 
scriptions of hadronization used in accelerator Monte 
Carlos. The Lund model is based on the 1980s picture 
of the link between QCD and string theory. The strings 



between quarks live in the 3+ld of QCD. The model as- 
sumes that a very long string forms between the quark 
anti-quark pair which then sequentially fragments with 
the fragments being assigned as various hadronic states. 
It's possible to morph our picture into the Lund one by 
assuming that the D7 brane in Figure 2 lies very close to 
the hard wall - the fifth dimension then plays little role 
and the string almost lies in the same plane as the end 
points move. How long the string grows and how many 
times it then breaks are not things we have computed so 
we can not dispute the Lund model. The striking differ- 
ence though between the models is that significant rho 
meson production occurs at the end points of the string 
moving in AdS, removing the need for repeated breakings 
of the initial string. Of course the differences in the AdS 
picture may be an artefact of a large N expansion and 
not relevant to true QCD, but equally N = 3 is believed 
to not be so far from N ^ oo. 

In theories close to the AT = 4 theory the rho mesons 
are special, in that they are associated with operators 
whose dimensions are protected from renormalization. 
This means they are present as supergravity modes in the 
DBI action of the D7 brane - other quark bound states 
would be represented by stringy states also tied to the 
D7 world volume. In QCD we would not expect a sepa- 
ration in character between these modes and all hadronic 
species should be produced at the end point governed by 
the same end point motion. 

The idea of separating, radiating string fragments is 
reminiscent of another reasonably successful model of 
hadronization. Thermal models [42 have been pro- 
posed that treat the event as separating fireballs radi- 
ating hadrons in thermal equilibrium. Some discussion 
of how that thermal spectrum can emerge from a gauge 
theory event with many final states even at zero temper- 
ature can be found in |43l |44] - it seems likely that the 
core idea is just that the energy of hadronization is freely 
available to all modes. A very toy model of holographic 
hadronization was presented, based on these ideas, in [45^ 
- there an equal lump of energy was dumped into the 5d 
holographic fields associated with each set of QCD bound 
states with a Gaussian profile in the radial direction. The 
model seemed to reproduce the data reasonably well too. 
j46] Here one can associate that lump of energy to the 
wave of gauge field emitted by the string end point as 
shown in Figures 9-11. 

It's therefore interesting that our holographic model 
seems to include aspects of both the Lund string model 
and thermal models of hadronization. We hope that in- 
sights from AdS will lead to phenomenologically more 
successful models of hadronization in QCD in the future. 



Acknowledgements: JF and ET would like to thank 
STFC for their studentship funding. NE is funded by an 
STFC rolling grant and is grateful to Marija Zamaklar 
and Kasper Peeters for discussions. KJ would like to 
thank Andreas Karch for many fruitful conversations. KJ 
was supported in part by the U.S. Department of Energy 
under Grant No. DE-FG02-96ER40956. 



17 

[9: 
[lo: 

[11 
[12: 
[13: 

[14 

[is: 
[16: 

[17 

[is: 

[19 

[20' 
[21 
[22: 
[23: 



Math. Phys. 2 (1998) 505 



J. M. Maldacena, Adv. Theor. Math. Ph ys. 2, 231 

_(1998) Int. J. Theor. Phys. 38, 1113 (1999) [arXiv ihep- 

th/9711200. 

[2] E. Witten, Adv. Theor. Math. Phys. 2 (1998) 253 

arXiv:hep-th/9802150 . 
[3] E. Witten, Adv. Theor. 

arXiv:hep-th/9803131 . 
[4] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, Phys. 

Lett. B 428 (1998) 105 arXiv:hep-th/9802109 . 
[5] M. Grana and J. Polchinski, Phys. Rev. D65 (2002) 

126005, [arXiv: hep-th/0106014 . 
[6] M. Bertolini, P. Di Vecchia, M. Frau, A. Lerda and 
R. Marotta, Nucl. Phys. B 621, 157 (2002) arXiv:hep-l 



th/0107057). 



I. Kirsch and D. Vaman, Phys. Rev. D72 (2005) 026007, 
[arXiv: hep-th/0505164 

A. Karch and E. Katz, JHEP 0206, 043 (2002) 
arXiv:hep-th/0205236 . 

M. Kruczenski, D. Mateos, R. C. Myers and D. J. Win- 
ters, JHEP 0307 049, 2003 ^arXiv :hep-th/ 0304032, . 
J. Babington, J. Erdmenger, N. J. Evans, Z. Guralnik 



and I. Kirsch, Phys. Rev. D 69 (2004) 066007 [arXivihep- 
th/0306018 . 

K. Ghoroku and M. Yahiro, Phys. Lett. B 604 (2004) 
235 ar Xiv:hep-th/0 408040 ■ 

J. Erd menger, N . Evans, L Kirsch and E. Threlfall, 
arXiv:071 1.4467 [hep-th]. 

S. D. Ellis, J. Huston, K. Hatakeyama and P. Loch, Prog. 
Part. Nucl. Phys 100 (2008) 484 arXiv:0712.2447| [hep- 
ph]]. 

L. F. Alday and J. Maldacena, JHEP 0711 (2007) 068 
arXiv:0710.1060 [hep-th]]. 

Y. Hatta, E. lancu and A. H. Mueller, JHEP 0805 (2008) 
037 arXiv:0803.2481 [hep-th]]. 

P. M. Chesler, K. Jensen and A. Karch, Phys. Rev. D 79 
(2009) 025021 arXiv:0804.3110 [hep-th]]. 
J. M. Maldacena, Phys. Rev. Lett. 80 (1998) 4859 
arXiv:hep-th/9803002 . 

S. J. Rey and J. T. Yee, Eur. Phys. J. C 22 (2001) 379 
arXiv:hep-th/9803001 . 

C. P. Herzog, A. Karch, P. Kovtun, C. Kozcaz 
and L. G. Yaffe, JHEP 0607 (2006) 013 arXivihep-j 
th/0605158 . _ 
_S. S. Gubser, Phys. Rev. D 74 (2006) 126005 |arXiv:hep-| 
th/0605182 . 

L. Girardello, M. Petrini, M. Porrati and A. Zaffaroni, 
JHEP 9812 (1998) 022 arXiv:hep-th/9810126 . 
K. Pilch and N. P. Warner, Nucl. Phys. B 594 (2001) 
209 arXiv:hep-th/0004063 . 

L. Girardello, M. Petrini, M. Porrati and A . Zaffaroni, 
Nucl. Phys. B 569 (2000) 451 |arXiv:hep-th/9909047, . 



[24] K. Pilch and N. R Warner, Adv. Theor. Math. Phys. 4 
(2002) 627 arXiv:hep-th/0006064 

[25] J. Polchinski and M. J. Strassler7arXiv:hep-th/0003136' 

[26] J. Babington, D. E. Crooks and N. J. Evans, Phys. Rev. 
D 67 (2003) 066007 |arXiv:hep-th/0210068 . 

[27] A. Kehagias and K. Sfetsos, Phys. Lett. B 454 (1999) 
270 arXiv:hep-th/9902125 



S. S. Gubser, arXiv:hep-th/9902155 



N. R. Const able and R. C. Myers, JHEP 9911 (1999) 
020 arXiv:hep-th/9905081 . 

J. Polchinski and M. J. Strassler, Phys. Rev. Lett. 88 
(2002) 031601 arXiv:hep-th/0109174 . 
J. Erhch, E. Katz, D. T. Son and M. A. Stephanov, Phys. 
Rev. Lett. 95 (2005) 261602 arXiv:hep-ph/0501128 . 
L. Da Rold and A. Pomarol, Nucl. Phys. B 721 (2005) 
79 arXiv:hep-ph/0501218 . 

L H. Brevik, K. Ghoroku and A. Nakamura, Int. J. Mod. 
Phys. D 15 (2006) 57 arXiv:hep-th/0505057 . 
J. Dai and J. Polchinski, Phys. Lett. B 220 (1989) 387. 
F. Bigazzi and A. L. Cotrone, JHEP 0611 (2006) 066 
arXiv:hep-th/0606059 . 

A. L. Cotrone, L. Martucci and W. Troost, Phys. Rev. 
Lett. 96 (2006) 141601 arXiv:hep-th/0511045 . 
K. Peeters, J. Sonnenschein and M. Zamaklar, JHEP 
0602 (2006) 009 arXiv:hep-th/0511044 . 
N. Evans and E. Threlfall, Phys. Rev. D 78 (2008) 105020 
arXiv:08 05.0956, [hep-th]]. 

K. Peeters, J. Plefka and M. Zamaklar, JHEP 0411 
(2004) 054 arXiv:hep-th/0410275 . 
[40] E. Poisson, Living Rev. Rel. 7 (2004) 6 [arXivigr- 
f qc/0306052. 

[41] B. Andersson, G. Gustafson, G. Ingelman and T. Sjos- 
trand, Phys. Rept. 97 (1983) 31. 
F. Becattini, Z. Phys. C 69 (1996) 485. 
J. Hormuzdiar, S. D. H. Hsu and G. Mahlon, Int. J. Mod. 
Phys. E 12 (2003) 649 arXiv:nucl-th/0001044 . 
Y. Hatta and T. Matsuo, Phys. Rev. Lett. 102 (2009) 
062001 |arXiv :0807.0098 [hep-ph]]. 

N. Eva ns and A. Tedder, Phys. Rev. Lett. 100 (2008) 
162003 |arXiv:071 1.0300 [hep-ph]]. 

Note to compare to data there is the need to impose 
an elaborate decay chain between the initially produced 
particles and the finally observed particles so the number 
of species for which there is data is much less than the 
total number of species in the initial yield which may 
hide many evils. 



